Combined method to extract spectroscopic information 
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Spectroscopic factors (SF) play an important role in nuclear physics and astrophysics. The tra- 
ditional method of extracting SF from direct transfer reactions suffers from serious ambiguities. 
We discuss a modified method which is based on including the asymptotic normalization coefficient 
(ANC) of the overlap functions into the transfer analysis. In the modified method the contribution 
of the external part of the reaction amplitude, typically dominant, is fixed and the SF is determined 
from fitting the internal part. We illustrate the modified method with (d,p) reactions on 208 Pb, 12 C, 
and 84 Se targets at different energies. The modified method allows one to extract the SF, which do 
not depend on the shape of the single-particle nucleon-target interaction, and has the potential of 
improving the reliability and accuracy of the structure information. This is specially important for 
nuclei on dripline, where not much is known. 

PACS numbers: 21.10.Jx, 24.10.-i, 24.50.+g, 25.40.Hs 



SF were introduced by the shell model formalism and 
are typically related to the shell occupancy of a state n 
in one nucleus relative to a state to in a nearby nucleus 
0. Today, phenomenological SF are extensively used 
in a variety of topics, from nuclear reactions to astro- 
physics or applied physics, yet the procedure for their 
extraction from the data has remained essentially the 
same for decades. For more than forty years since the 
dawn of nuclear physics, direct transfer reactions, such 
as (d,p), (d,t), ( 3 He, dY, He, a), have been the central 
tool to determine SFs H,y,Q. Extracting SFs with good 
precision from data is very important to test the valid- 
ity of today's many body theories. For conventional nu- 
clei there are many experiments available providing SFs, 
which are often lower than those predicted by shell model 
0. Electron-induced knockout or electron scattering is 
supposed to provide a better accuracy in extracting SFs 
than transfer However, for exotic nuclei near or 

on the driplines, transfer reactions are a unique tool and, 
hence, can have a large impact in the programs of the 
new generation rare isotope laboratories. Given the ex- 
perimental difficulties faced with measurements on the 
driplines, it is crucial to have a reliable method for ana- 
lyzing and extracting useful information from each single 
data set. 

Usually, transfer angular distributions are analyzed 
within the framework of the distorted- wave Born approx- 
imation (DWBA). The SF determined by normalizing 
the calculated DWBA differential cross section to the 
experimental one (e.g. 0, S 0) is compared with the 
SF predicted by shell model. Even when error bars in 
the experimental cross section are low, the uncertainty of 
the extracted SF resulting from the normalization of the 
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DWBA cross section is often large, regardless of whether 
it agrees with the shell model prediction. The reasons for 
this inaccuracy are typically: i) optical potentials ambi- 
guity, ii) the inadequacy of the DWBA reaction theory, 
or iii) the dependence on the single-particle potential pa- 
rameters. The first point has been object of a recent 
systematic study ^(J- The second point needs to be ad- 
dressed case by case, and examples of improved reac- 
tion models are the coupled channel Born approximation 
(e.g. llj) or the continuum discretized coupled chan- 
nel method (e.g. |12|). This work will critically review 
the standard procedure of extracting SFs from trans- 
fer reactions focusing on the third point; the modified 
method eliminates the dependence of the extracted SFs 
on the single-particle potentials, the main advantage of 
the method. 

We will address a modified approach to spectroscopy 
from transfer reaction which includes the asymptotic nor- 
malization coefficient (ANC) in the analysis 0. For sim- 
plicity, in the following formulation, we consider A(d,p)B 
reaction and disregard spins (naturally these are included 
in the applications). The DWBA amplitude for this re- 
action is given by: 

M=<ip { f ] lAn\^V\ip pn 4 +) >, (1) 

where AV = V pn + VpA — U P B is the transition operator in 
the post-form, Vij is the interaction potential between i 
and j, U P B is the optical potential in the final-state. The 
distorted waves in the initial and final states are 
and tp j, ^ , ip pn is the deuteron bound-state wave function 

and (r) is the overlap function of the bound-states of 
nuclei B and A which depends on r, the radius-vector 
connecting the center of mass of A with n. The overlap 
function is not an eigenfunction of an Hcrmitian Hamil- 
tonian and is not normalized to unity [l3|. The square 
norm of the overlap function gives a model-independent 
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definition of the SF: 



S = N </M„> 



(2) 



Here, N is the antisymmetrization factor in the isospin 
formalism (N will be included in the overlap function 
from now on). 

The leading asymptotic term of the radial overlap func- 
tion (for B = A + n) is 



llniij) ( r ) ^ °ij i*hi(iKr), 



(3) 



where hi(inr) is the spherical Bessel function, n = 
\/2 [i An £An, £ An is the binding energy for B — > A + n, 
and HAn is the reduced mass of A and n. Similarly, 
the asymptotics of the neutron single-particle wave func- 

r>R 

tion is <fiAn(n r ij) ( r ) ~ b nr ij i nh[(i nr), where n r is the 
principle quantum number. The asymptotic behaviour is 
valid beyond R, the channel radius. It is clear that, in the 
asymptotic region, the overlap function is proportional 
to the single particle wave function. The normalization 
Cij introduced in Eq. J2J is the ANC which relates to the 
single-particle ANC (SPANC) 6„ r y by Cy = K nrlj 6 nr y, 
where K Ur ij is an asymptotic proportionality coefficient. 
It is standard practice to assume that the proportional- 
ity between the overlap function and the single particle 
function extends to all r values 



An(y)( r ) = Kn r lj ¥>An(n r lj){r)- 
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Since fAntnAj) ( r ) ^ s normalized to unity, this approx- 
imation (Eq. 0} implies that Sy = K^ u- Wc 
have to emphasize, however, that the overlap func- 
tion in the interior is nontrivial and may well dif- 
fer from the single particle wavefunction. Approx- 
imating the radial dependence of the overlap func- 
tion as described above leads to the DWBA ampli- 
tude M = K„ r ij < ipf<PAn(.nrlj)\& v \'PpnW >■ Nor- 
malizing the calculated DWBA cross section, 
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I < <PAn(n r lj)\AV\(p pn l{j. 
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to the experimental data provides the phenomenological 
SF Sij = Kn r lj- Assuming that Eq.QJ is valid for all r, 
we can infer from Eq.J2J that the main contribution to 
the norm of the overlap function comes from the nuclear 
interior. 

In order to make the dependence on the SPANC more 
explicit, we split the reaction amplitude into an interior 
part and an exterior part: 



M 



K nr ij Mint [6] + Knrlj bn r lj M ext , 



(6) 



where the internal part of the matrix element 

M in t[6„ r y] =< ^f ) ip An (n r lj)\^y\ l Ppni'i >r<R 
depends on b nr ij through the bound state wave- 
function tpAn(n r lj) i while the external part 
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ext =< ipf ■* inhi{i Kr)\AV\ippn ipl +> > r >R does 
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not depend on 6 nr y. Here, R is the channel radius 
taken so that for r > R the overlap function can be 
approximated by its asymptotic form Eq.© (R is 
only used to illustrate the method as in the end this 
separation is not required). The contribution from the 
nuclear exterior is fixed by the ANC, whereas the SF 
determines the normalization of the internal part of 
the radial matrix element. Since transfer reactions are 
dominantly peripheral, SFs can only be extracted from 
transfer reactions due to a small contribution from the 
nuclear interior. We now introduce the ANC into the 
DWBA cross section: 
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Introducing Eq.© into Eq.(© and dividing by Cfj, we 
arrive at a function R DW (b) 



R DW {bn rl] ) = | 
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Note that the single-particle ANC 6„ r y itself is a function 
of the geometrical parameters of the bound state n — A 
nuclear potential (ro, a) which are, a priori, not known. If 
the ANC and the cross section for the (d,p) reaction have 
been measured, the experimental counterpart of R DW , 
R ex P = d j n / Cfj can be experimentally fixed. Then, 



(4) imposing the equality 



R exp = R 
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(bn r lj), 
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will provide the correct 6„y and consequently the SF 

sij = cybi lr 

At this stage, a few points should be made clear. First 
of all, for specific optical potentials, Eq. (J5J depends on 
two independent parameters, Sy and 6„y . In the stan- 
dard approach, to evaluate this cross section, the second 
parameter is fixed by arbitrarily choosing the bound state 
n — A potential geometry. Thus, the extracted product 
Sij b\ y does not coincide necessarily with the correct 
ANC. Since the ANC determines the normalization of 
the external part of the DWBA amplitude, in the stan- 
dard approach the SF is determined by an unrealistic 
variation of the external contribution. In the modified 
method here discussed, since the contribution of the ex- 
ternal part is fixed through the correct ANC, the whole 
DWBA procedure loses this artificial degree of freedom. 

Secondly, if the reaction is peripheral, i. e. the first 
term in Eq. JHJ is negligible, one can determine the ANC. 
So, the modified approach makes use of two experiments: 
the first to fix the ANC, the second to determine the SF 
consistent with that ANC. In present experiments and 
with the new generation of rare isotope facilities, ANCs 
can be determined with 5% accuracy. Since the deter- 
mination of the SF comes from the internal region, the 
second experiment needs to be performed at a beam en- 
ergy for which the contribution from the interior is signif- 
icant. The higher the contribution of the internal region, 
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FIG. 1: Cross section for 208 Pb(d,p) 209 Pb(g.s.) at 22 MeV 
and the dependence on the single particle parameters: ex- 
perimental value (solid line), experimental error bar (shaded 
area) and the DWBA prediction (dot-dashed). 

the stronger the dependence on b nr ij in R DW (b nr ij) and 
the smaller the uncertainty of the extracted SF, although 
a balance needs to be found since large interior contribu- 
tions may not be well describe by DWBA. The DWBA 
differential cross section near the main peak of the angu- 
lar distribution and, correspondingly, R DW (b nr ij) are the 
functionals of the single-particle ANC b nr ij. One given 
bn r ij can be produced by an infinite number of single- 
particle potentials, local and non-local. However, the 
dependence of da DW /dfi or R DW (b„ r ij) on the shape 
of the single-particle potential is minor. Hence, the ex- 
tracted SF in the modified method does not depend on 
the single-particle potential. We illustrate the method 
presenting three different applications: i) 209 Pb, ii) 13 C 
and iii) 85 Se. We will drop the subscripts on b for sim- 
plicity. 

Let us consider the reaction 208 Pb(d,p) 209 Pb from [l^ . 
Although the ANC for < 209 Pb| 208 Pb> is not published, 
it can be determined from the sub-Coulomb reaction [LH ] 
208 Pb( 13 C, 12 C) 209 Pb as the other vertex < 13 C\ 12 C > is 
well known [l5j. Sub-Coulomb reactions are extremely 
peripheral and insensitive to details of the optical po- 
tentials. For this reason they present an excellent probe 
for extracting the ANC accurately. From [l4| we obtain 
an ANC C 2 g/2 = 2.15(0.16)fm- 1 for 209 Pb. Then using 

208 Pb(d,p) 209 Pb data at E d = 22 MeV we obtain 
R ex P = 2.46(0.31) fm mb/srad, where the error bar is 
calculated based on both, the ANC and the cross sec- 
tions errors, taken as independent. The experimental 
data in [12j has 1% accuracy but is taken only down to 
6cm = 35° whereas the peak of the DWBA distribution 
is at 8 cm — 25°. We extrapolate the data based on the 
shape predicted by DWBA and include a 10% error in 
the cross section to account for this difference. Measure- 
ments at 25° could improve the error bar in R ex P consid- 
erably. We next perform a series of finite range DWBA 
calculations for 208 Pb(d,p) 209 Pb (E d = 22 MeV), using 
the optical potentials from 0. The adiabatic prescrip- 
tion was used to take into account deuteron breakup 
which is important for this reaction. The Reid-soft-core 
potential was used for the deuteron wavefunction, as well 
as in all other examples. For illustration purposes, we 




FIG. 2: Cross section for 12 C(d,p) 13 C(g.s.) at 51 MeV and 
the dependence on the single particle parameters: experimen- 
tal value (solid line) , experimental error bar (shaded area) and 
the DWBA prediction (dot-dashed). 

use a Woods Saxon well to generate the 208 Pb + n single- 
particle wavefunctions and obtain a range of SPANCs b 
by varying the single particle parameters (tq , a) and ad- 
justing the depth to reproduce the correct binding for the 
2.t/ Q/9 in each case. We use the same s.o. strength as that 
in 14] although the s.o. strength does not affect the final 
result. 

The results of our calculations R DW (dot-dashed line) 
and the experimental value R ex P (solid line and shaded 
area) are presented, as a function of b, in Fig. ^ From 
R ex P one finds b = 1.82 fm" 1 / 2 and S=0.74. It is worth 
noting that in the standard approach typical parame- 
ters (r ,a) = (1.2,0.6) fm, produce b = 1.34 fm" 1 / 2 . 
The direct comparison of the DWBA cross section using 
(ro,a) = (1.2,0.6) fm, with the data, give S=0.866 and 
consequently, C 2 = 1.56 fm -1 , beyond the experimental 
range. As pointed out before, in the standard approach 
the SF is determined at the cost of an artificial ANC. 

The beam energy of 22 MeV is above the Coulomb 
barrier, thus the reaction is not peripheral. This can be 
seen in Fig. ^through the slope of the dot-dashed curve. 
In fact for this particular energy, the interior contribution 
is around 10%. The uncertainty in b 6 [1.1,3.1] fm -1 / 2 
propagates into a large uncertainty in S € [0.3, 2.2]. This 
is due to the fact that the contribution from the interior 
at this energy is still small. The scaling factor relating 
the uncertainty of S with that of b is (|M> \/\M < |) 2 . The 
smaller the contribution from the interior, the smaller 
the accuracy with which the SF can be determined. 

Also in Fig. ^ we show the results for R DW correspond- 
ing to the calculation at Ed — 5 MeV (dashed line) . This 
is to illustrate that, at sub-Coulomb energies, the reac- 
tion becomes completely peripheral and the dependence 
on b disappears. Measurements at these energies could 
provide C| with accuracy < 5%. In addition, measure- 
ments at higher energy (> 30 MeV) would increase the 
slope of R DW (b) and decrease further the error on the 
extracted SF. 

Another standard case is the 12 C(d,p) 13 C reaction, 
for which many data sets are conveniently compiled in 
a recent publication 0. We studied three cases (8.9 
MeV, 30 MeV and 51 MeV), using the same JLM op- 
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tical potentials as 10]. We perform a series of finite 
range DWBA calculations varying the lpi/2 12 C-n sin- 
gle particle parameters, in order to obtain R DW (b) as 
described before. Results for the less peripheral case (51 
MeV) are plotted in Fig. H (dotted-dashed line). We 
take the data from and the ANC from [l^ , to obtain 
Rexp = a{2^ = 2.92(0.35) fm mb/srad. An S = 0.66 

1,1/2 

(shell model) would require b = 1.89 which is contained 
in our results. However, such a conclusion is misleading. 
Fig. |3 shows that even for this relatively large energy 
the dependence of R DW on b is weak. Consequently, it 
is not possible to extract a SF. 

It was pointed out in [T^j that the deuteron breakup 
is important for this reaction and should be taken into 
account. To emphasize this fact, we compare our re- 
sults using the adiabatic deuteron potential 0] from [ljj 
(dot-dashed line in Fig. [5J with those obtained using an 
optical potential fitted to the deuteron elastic scattering 
(dotted line in Fig. EJ. The disagreement is very large. 
Interestingly, the method here described is also able to 
detect inadequate optical potential parameterizations. 

Oak-Ridge has developed a program to measure a se- 
ries of inverse kinematics (d,p) reactions for nuclei on the 
neutron dripline |18| . As one of the nuclei in the program 
is 85 Se, we have performed exploratory calculations for 
84 Se(d,p) 85 Se. We take global parameterizations for the 
optical potentials and perform a series of calcula- 
tions varying the single particle parameters. We com- 
pare the dependence of R DW on b for a range of energies 
E d = 4 - 100 MeV. We verify that, expectedly, the de- 
pendence on b increases with beam energy. We find that 
Oak-Ridge energies (10 MeV/A) are adequate to deter- 
mine ANCs but not SFs. However, a facility that allows 
for the production of 84 Se at E > 25 MeV/A (such as 
NSCL-MSU, GANIL or RIKEN) could provide accurate 
spectroscopic information. 

In conclusion, we have presented an alternative method 
to extract SFs, taking into account the sensitivity of the 
transfer data to the interior part of the overlap function 
and combining that information with the ANC. Transfer 
data can only become useful within this method if it has 
a significant contribution from the interior, and is well 



described through a one-step DWBA formalism. The 
balance between these two conditions is not a trivial one. 
By reducing the error bars in both the measured transfer 
cross section and the ANC, this prescription determines 
the single particle asymptotics and from it, a SF with 
reduced uncertainty. The ANC needs to be determined 
independently; it can be pinned down accurately with the 
same transfer reaction at sub-Coulomb energies or using 
heavy-ion induced reactions, both safely peripheral. Note 
that uncertainties due to optical potentials and higher 
order effects need to be assessed independently, as this 
work focuses on the single particle parameter uncertain- 
ties only. 

The method here presented has the potential of re- 
ducing the uncertainty in the overlap function consider- 
ably. However it still assumes that the interior part has 
a Woods Saxon single particle wavefunction shape. This 
has been corroborated by recent Green's Function Monte 
Carlo calculations on light nuclei [19( . Even if there were 
non-localities of the single particle potential this would 
affect mostly the deep interior and thus would not be 
visible in the transfer reactions. 

Results for (d,p) on 208 Pb were used to illustrate the 
method. We discussed previous analyzes of (d,p) reac- 
tions on 12 C, and showed the limitations. We have also 
demonstrated that this method can rule out inadequate 
choices of optical potentials. Considering specific future 
experiments, we have performed exploratory calculations 
for (d,p) on 84 Se. This method will become useful for a 
broad variety of transfer experiments in the field of rare 
isotopes. The same method can equally be used for trans- 
fer to excited states. These same ideas can be extended 
to other reactions, in particular breakup reactions which 
also have an impact on Astrophysics. Finally, it would 
be helpful if the state-of-the-art reaction codes would in- 
corporate the formalism discussed. 
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